Precision Test of Mass Ratio Variations with Lattice-Confined Ultracold Molecules 
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We propose a precision measurement of time variations of the proton-electron mass ratio using 
ultracold molecules in an optical lattice. Vibrational energy intervals are sensitive to changes of 
the mass ratio. In contrast to measurements that use hyperfine-interval-based atomic clocks, the 
scheme discussed here is model-independent and does not require separation of time variations 
of different physical constants. The possibility of applying the zero-differential-Stark-shift optical 
lattice technique is explored to measure vibrational transitions at high accuracy. 



Ultracold molecules open new opportunities for preci- 
sion measurements of possible variations of fundamen- 
tal physical constants. The test of time variation of the 
proton-electron mass ratio A/i//i (/i = mp/rrie, where 
•nip and me are the proton and electron masses) is par- 
ticularly suitable, since molecules arc bound by electronic 
interactions while ro- vibrations are dominated by nuclear 
dynamics. Recent proposals to search for A/i//i include 
detecting changes in the atomic scattering length near 
a Feshbach resonance [l] or using near-degeneracies of 
molecular vibrational levels from two different electronic 
potentials to probe microwave frequency shifts arising 
from A/i/^ 0, 0]- We propose a two-color (Raman) op- 
tical approach to directly determine vibrational energy 
spacings within a single electronic potential of ultracold 
dimers in an engineered optical lattice. The measure- 
ment relies on the cumulative effect of A^//^ on the ex- 
cited vibrational levels, and utilizes the entire molecular 
potential depth to enhance precision by choosing two vi- 
brational levels with maximally different sensitivity. 

For a given physical system, we can write a propor- 
tionality relation between A/i//i and the corresponding 
fractional change in transition frequency l^v/v as 



A^ 



(1) 



Only dimensionless quantities appear in Eq. ([T]), to avoid 
introducing a time dependence of the Cs hyperfine fre- 
quency, for example. In the case discussed here, n is of or- 
der one and has a small potential-dependent uncertainty. 
The fractional uncertainty 5^i/ ^. of the measurement of 
(A/x ± Sfj,)/ fj, must be minimized. For a given frequency 
measurement uncertainty Sv, from Eq. JT]) we obtain 
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The last step in Eq. ([2|) is motivated by the assumption 
that experimental limitations constrain 5v rather than 
5v/v. Eq. ([2]) indicates that the quantity (di^/dlnfi) 
must be maximized. In other words, we search for the 
energy gap with the maximum absolute frequency shift 
arising from a given fractional mass ratio change. While 



a microwave measurement [3, Isj could have a smaller Su 
than the optical frequency-comb-based Raman approach, 
the latter maximizes sensitivity through the cumulative 
effect of the entire molecular potential depth. 

Heteronuclear dimers may be advantageous for the 
present proposal due to deeper electronic ground state 
potentials by a factor of ^^2-5 {dv/ dlnfi is proportional 
to the same factor). On the other hand, homonuclear 
dimers can lead to higher precision since their radiatively 
long-lived vibrational levels in the electronic ground state 
are insensitive to blackbody radiation. Molecules based 
on even isotopes of alkaline-earth- type atoms (e.g. Sr, 
Ca, Yb) enjoy the lack of hyperfine and magnetic struc- 
ture in the electronic ground state, simplifying the prepa- 
ration of the system and reducing systematic shifts. 

As in the work on a Sr-atom optical lattice clock Ij] and 
narrow-line photoassociation (PA) spectroscopy [5|, we 
propose to use ultracold Sr dimers confined in an optical 
lattice for the A/i//i experiment. Vibrationally excited 
Sr2 in the electronic ground state is produced via PA. 
Raman spectroscopy aided by a femtosecond optical fre- 
quency comb will be used to interrogate the energy spac- 
ings between deeply bound vibrational levels and those 
closer to the dissociation limit. The Franck-Condon fac- 
tors (FCFs) between the electronic ground state X poten- 
tial (dissociating to ^5*0 S'o) and the excited 0+ poten- 
tial (dissociating to ^Sq -\-^ Pi, with ungerade symmetry 
and the atomic angular momentum projection onto the 
molecular axis = 0) are sufficiently favorable to enable 
Raman transitions between two vibration levels in the X 
state of Sr2 via '^Pi, as both potentials are dominated 
by van der Waals interactions. Since the X potential is 
30 THz deep [1, 0, Hi, and the relative stabifity of the 
Raman lasers can be maintained to better than 0.1 Hz 
via the comb \^ , we expect a precision of better than ^ 
(O.I Hz)/(10 THz) = 10'^^ in the test of Afi/fi. 

Other tests based on atomic frequency metrology con- 
strain A/i//i to ^ 6 X 10^^^/year and the recent 
evaluation of astronomical NH3 spectra limits A/i//i to 



3 X 10 ^^/year [llj. The atom-based tests rely on 



theoretical interpretations such as the Schmidt model, 
since electronic and fine structure transitions do not di- 
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FIG. 1: (a) The model potential curve for the Sr2 ground state 
(solid line) , and the Morse potential fitted to three parameters 
(dashed line), (b) Vibrational energy sensitivities to Afi/fi, 
as a function of the vibrational number v. 




FIG. 2: The proposed Scheme I for precision Raman spec- 
troscopy of Sr2 ground state vibrational level spacings. A 
two-color photoassociation pulse prepares molecules in the 
V — —3 vibrational level. Subsequently, a Raman pulse cou- 
ples V — —3 and v = 27, via v' ~ 40 of 0^. 



rectly depend on and hyperfine transitions simultane- 
ously depend on /i, a (the fine structure constant), and 
i?oo (the Rydberg constant). The NH3 result is based 
on molecular lines and is therefore less model-dependent, 
but relies solely on cosmological observations. In addi- 
tion, it disagrees with the H2-based result relevant to the 
same cosmological age (10"'^'' years) that indicates non- 
zero A/i//i at the 10~^^/year level The molecular 
system proposed here provides a direct test of present- 
day variations with a competitive precision and a weak 
dependence on theoretical modeling [13]. 

In this work we select vibrational levels of the ground 
state potential that have the largest and smallest sensi- 
tivities to A/z//x. To model the Sr2 ground state poten- 
tial V{r), we combine the experimental RKR potential 
01 with its long-range dispersion form —ce/r^ — cs/r^, 
where cg = 3100 Eha% jlil, cg = 1.9 x 10^ Eha% is de- 
terminted by smoothly connecting to the RKR potential, 
and r is the internuclear separation in units of {Eh — 
4.36 X 10"^^ J and as = 0.0529 nm). The model poten- 
tial has 61 vibrational levels, depth d = 4.7 x 10^^ Eh, 
the minimum at the internuclear separation tq = 8.9 as, 
and a scattering length of ^ 8 ob JJ, 18]. The Morse po- 
tential can be quantized analytically and is a convenient 
approximation to the ground state molecular potential, 

VMir)=d{l-e--^---«^)^ -d, (3) 

where a ~ 0.7 a^^ for Sr2. The Morse energy levels are 

en = 2eoin+^)-elin+^f/d-d, (4) 

where eo is approximately the zero-point energy. Note 
that the Morse spectrum is valid only if e„ — Cn-i > 0, 
or n < d/eo, which means that Vm('') has about N ~ 



40 bound levels. Figure [T] (a) compares the Sr2 model 
potential and its Morse approximation. 

Since eo oc fJ,^^^^, the logarithmic derivative of the 
expression for the nth vibrational energy level is 

^^-e.(.. + i) + .^(n + i)V<i. (5) 

Equation ([5]) is the energy level sensitivity to a fractional 
mass ratio change, with the maximum absolute sensi- 
tivity for TT-niax — N/2, and lowest sensitivity near the 
bottom and the top of the potential well, as expected 
given a fixed potential depth. The sensitivities were also 
determined for our Sr2 model by calculating vibrational 
energies for slightly different atomic masses. Both the 
Morse approximation and more realistic calculation point 
to 25 < w ^ 28 as the most sensitive. Figure [1] (b) shows 
the level sensitivities to A/i//z. We choose to work with 
V — 27, and the Afi/ 11 measurement is optimized if either 
a weakly bound or the deepest vibrational level is chosen 
as the reference level for Raman spectroscopy. 

Based on considerations of sensitivity as well as molec- 
ular transition strengths, we propose two schemes for the 
measurement of Afi/ fj,. While they yield the same infor- 
mation, their combination provides a consistency check 
on the experimental method. The first scheme is more 
straightforward as it relies on one of the least-bound vi- 
brational levels for the Raman transition. The second 
scheme involves an additional Raman step to drive the 
weakly bound molecules into a deeper vibrational level. 
Below, the rotational angular momenta are J = for 
X and J' = 1 for 0+ and 1„. The transition strengths 
are obtained from relativistic configuration interaction ab 
initio calculations [H, [13] adjusted to agree with mea- 
surements of weakly bound 0+ vibrational levels Q . 

Scheme I is illustrated in Fig. [21 and measures the 
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FIG. 3: Calculated Fran ck- Condon factors, or transition 
dipole moments squared, between vibrational levels of the X 
and electronic states as a function of the binding energy 
for vibrational levels of 0^ (J = and J' — 1). As relevant 
to the text, only the FCFs of n = 0, 27, and —3 are shown. 
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FIG. 4: Magic frequencies for the optical lattice near 910 
nm (10990 cm~^). The detunings from resonance are ~ 1.5 
and —3 GHz for Schemes I and II, respectively, and must be 
confirmed experimentally. 



energy difference between the weakly bound v = —3 
(negative vibrational numbers imply counting from the 
top of the potential with the least-bound level being —1) 
and V = 27, the latter being most sensitive to mass ratio 
changes. Step la is two-color PA into v = —3 via v' = —6 
using 689 nm light, with the two colors detuned by 
about 10 GHz (the primes refer to vibrational levels of 
the 0„ excited electronic potential) . Step lb is three-level 
Raman spectroscopy, v = —3 ^ -y' ~ 40 ^ v = 27. Fig- 
ure [3] shows the FCFs, that include the transition dipole 
moments and are defined as \{v, J\er\v' , J')P, for v = 0, 

V = 27, and v = —3 to any vibrational level of 0+. For 
?; = — 3 the FCFs approach lO"** (eos)^ for a number of 
0^ vibrational levels with binding energies smaller than 
1000 cm~^, where e is the electron charge. For v — 27 
the maximum FCFs are about hundredfold larger, and 
quickly decrease for binding energies smaller than 400 
cm"-'^. This suggests using 0^ intermediate levels with 
binding energies around 400 cm~^ to balance the Raman 
transition strengths, such that the laser wavelengths are 
in the 700-750 nm range. The resulting FCFs of ~ lO"'' 
(ea^)^ imply that for a 1 GHz detuning from the inter- 
mediate level v' ^ 40 and laser intensities of 2 W/cm^, 
the two-photon Rabi rate is '^-^ 27r x 10 Hz. 

Scheme H that probes deeper levels is possibly less sen- 
sitive to collisional relaxation. Its disadvantage is an ex- 
tra step of two-photon population transfer. Step Ila is 
the same as la. Step lib is Raman population transfer 

V = —3 —> v' ^ 40 V — 27, analogous to Step lb. 
Finally, step He is w = 27 ^ w' ~ 10 w = Raman 
spectroscopy, with the wavelengths in the 700-800 nm 
range. As shown in Fig. [3l the FCFs between v = and 
the vibrational levels of 0^ are very large near the bottom 
of Oj and decrease rapidly for levels with binding ener- 
gies smaller than 1500 cm^^. This suggests the choice of 
v' ^ 10 for the intermediate level as it balances the FCFs 
within the Raman transition to ~ 10~^ (eas)^- For these 



transition strengths, and using 50 MHz detunings from 
v' ~ 10 and Raman laser intensities of 10 W/cuP, the 
Rabi rate for step He is ^ 27r x 10 Hz. 

The experiment critically depends on the control over 
systematic effects and the ability to perform Raman spec- 
troscopy on a large number of molecules for a good sig- 
nal to noise ratio. Using an optical lattice to trap the 
ultracold molecules Q is beneficial both for attaining 
high densities on the order of 10^^ /cm^, and for control- 
ling systematic shifts. The zero-differential-Stark-shift 
(or magic frequency) technique allows precise and ac- 
curate neutral-atom clocks [1, [3, It relies on the 
crossing of dynamic polarizabilities of the two probed 
states at a certain lattice frequency. Such a lattice en- 
sures a vanishing differential Stark shift and a suppres- 
sion of inhomogenous Stark broadening. For PA at the 
^5*0 -|-'^ Pi dissociation limit, the polarization-dependent 
magic wavelength is near 914 nm (10950 cm^^). 

Analogously, we can search for a zero-differential- 
Stark-shift lattice frequency for the proposed pairs of Sr2 
vibrational levels of the X potential. The Stark shifts 
of vibrational levels are proportional to dynamic polariz- 
abilities [ll[ll, and are independent of the light polar- 
ization for J = 0. The polarizability of Sr2 in the vicin- 
ity of 914 nm slowly decreases with vibrational quantum 
number of X. However, tuning the lattice frequency to 
near resonance with a vibrational transition from X to 
1„ makes it possible to match the polarizabilities of two 
vibrational levels in X. As shown in Fig. lU our calcula- 
tions indeed reveal lattice frequency values in the vicinity 
of 11000 cm~^ for which the polarizabilities of u = 27 and 
V — —3, as well as those of w = and v = 27, are equal. 
The precise location of the resonance will shift as more 
experimental input becomes available. 

A possible disadvantage of working near these reso- 
nances is enhanced scattering of lattice light. We esti- 
mate the scattering rate from the resonance strength as 
well as from the lattice detuning needed to achieve the 
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zero-differential-Stark-shift condition. From our model, 
these detunings are ~ 1.5 GHz for Scheme I and ~ —3 
GHz for Scheme H. Calculations show that the total 
spontaneous decay rate of 1„ levels with 4000 cm~^ bind- 
ing energy (as in Fig. S]) is I-k^ = 27r x 9 kHz. The ef- 
fective photon scattering rate is given by F^ = 27r7s = 
(7rs7) 7(25/7)^, where 6 is the detuning from the reso- 
nance. The measure of the lattice intensity / is s, such 
that s = ///sat, with the effective saturation intensity 
for a single vibrational channel /sat — ('''c?i^7^)/(4/Fc), 
where 27r?i is the Planck constant and /fc is the FCF 
(/fc ~ 5 X 10"^ (eas)^ for the resonance in Fig. [4]). The 
estimated scattering rates are thus ~ 4/s for Scheme I 
and ~ 1/s for Scheme H, for / = 10 kW/cm^. This in- 
tensity supports a trap a few /iK deep. Note that the 
scattering of lattice photons due to Oj vibrational levels, 
as well as non-resonant In levels, is strongly suppressed. 

Further, we estimate the incoherent scattering rate of 
Raman lasers by the intermediate 0+ levels. For the v — 
—3 — > 40 transition, for example, this scattering 

rate is < G.l/s for a 1 GHz detuning and 2 W/cm'^ laser 
intensity. For the u = 27 -> u' ~ 10 transition, it is 
~ 1/s for a 50 MHz detuning and 10 W/cm'^ intensity. 
These estimates are conservative since the excited state 
population is expected to be suppressed when the Raman 
condition is fulfilled. Thus the scattering rates of the 
lattice and spectroscopy lasers will not limit the ^ 10 Hz 
power-broadened linewidth of the two-photon transition. 

The Stark shift of v due to one of the Raman lasers is 
A = (s7^)/(8(5) ~ ^si&ll) in the large detuning limit. 
We estimate the near-resonant contributions to Stark 
shifts to be ~ 50 Hz for both schemes for each Ra- 
man laser. Moreover, the shifts of the two vibrational 
levels within a Raman pair have the same sign which 
leads to cancellation with the proper power balance of 
the Raman beams. In addition to the near-resonant shift, 
there is a background shift due to other molecular levels 
(which can be partially compensated by slightly shifting 
the power balance in the Raman beams). In the 700-800 
nm wavelength range and for Raman laser intensities of 
^ 10 W/cm^, the background Stark shifts are ~ 100 Hz. 
Hence, in the worst case the Raman beam intensity must 
be controlled to < 1% for a sub-Hz measurement. 

Other major systematic effects are magnetic field fluc- 
tuations, and Sr density variations in the lattice. Our 
past work on the ^^Sr atomic clock demonstrates 
that these effects can be controlled to below 1 Hz. Impor- 
tantly, the absence of magnetic structure in the ground 
electronic state of Sr2 with J — Q should significantly 
reduce any magnetic shifts. 

To experimentally search for the X vibrational levels, 
we will proceed in two steps. The first is to locate the 
ti = —3 level of X as follows. PA into v' — —6 leads 
to preferential decay into v — —3 Tuning a laser to 
V = —3 — !■ w' = —6 will convert the v — —3 molecules 



back to atom pairs. Hence, observation of a reduced atom 
trap loss associated with PA should allow locating the 
V — —3 level. Once the v = —3 level is identified, similar 
methods can be used to find the v' ~ 40 level of 0+ and 
the V = 27 level oiX, as well as the other proposed levels, 
in a stepwise manner. For example, if the molecules are 
transferred from v = —3 to v' 40, a reduced fraction 
of Sr2 will be converted back to atoms in the procedure 
outlined above. The lifetime of the v = —3 molecules 
should be long due to the lack of radiative decay channels. 
The coUisional relaxation rates are not well known and 
must be determined experimentally. 

In conclusion, our analysis of Sr2 dimers shows that ul- 
tracold non-polar molecules in a zero-differential-Stark- 
shift optical lattice is an excellent system for measuring 
time variations of mass ratios. This system provides a 
model-independent test that is based on different physics 
than atomic clocks. We expect a sub-Hz frequency mea- 
surement for a ~ 10^^'^/year test of A^/^, with future 
improvements in precision by at least a factor of ten. 
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